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ASYMPTOTIC  PHOBIAS  IB  TRAMBOHIC  PLOW 

Julian  0.  Cola 
Department  at  Mathematics 
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1.  Introduction 


This  lecture  Is  devoted  to  « .tuples  of  the  use 
of  asymptotic  matching  In  transonic  flow  theory. 

The  discussion  is  based  on  transonic  small -distur¬ 
bance  theory  (ISC)  whose  usefulness  Is  already  well 
understood.  By  the  use  of  asymptotic  expansions 
valid  in  different  regions  greater  utility  can  be 
given  to  several  exact  ISC  solutions,  such  as 
those  of  Mach  number  one,  or  strictly  two  dimen¬ 
sional  flow.  These  exact  solutions  can  be  used  as 
the  principal  part  of  a  flow- field  and  corrections 
can  be  calculated. 

In  the  applications  of  the  method,  typically, 
different  asymptotic  expansions  based  on  different 
limit  processes  are  constructed.  These  expansions 
are  valid  in  different  regions  of  space  but  can  he 
matched  asymptotically  by  a  class  of  intermediate 
limits  valid  in  an  overlap  domain  [1].  Symboli¬ 
cally  we  might  have 

f(x;e)  »  QQ(e)  f*(x*)  +  a^e )  f*(x*)  c<1«  aQ 

+  •  •  •  x*  -  x*(x ;  e) 

*  &0(e)  f0(x)  +  P^UjfjU) 

+  •••  x  =  x(x,e) 

The  associated  limits  have  e  -  0,  x*  fixed  as 
t-0  x  fixed  respectively.  x*,j  are  definite 
functions  of  (x.t).  Generally  speaking  each  ex¬ 
pansion  is  valid  in  a  different  region,  for  example 
x*  >  x^,  x  <  x„ .  Where  these  regions  are  adjacent 
then  the  validity  of  the  expansions  can  usually  be 
extended  to  an  overlap  domain  In  which  the  two 
expansions  must  agree  to  a  certain  order.  A  class 
of  limits  intermediate  to  the  two  above  is  expressed 
by  x,  «  x^  (x,e)  fixed  as  6-0  such  that  under 

this  limit  as  £  -  0,  x*  «  x*  (x(x. )  ,6;  -  0, 
x  »  x(x(x  i)  ,c)  -  ».  Utader  the  class  of  intermed¬ 
iate  limits  the  two  expansions  read  the  same  in 
x  ;  thus  information  about  constants  of  integra¬ 
tion  and  the  form  of  the  expansions  can  be  found. 
Various  applications  of  this  idea  in  a  broader  con¬ 
text  appear  in  the  lecture. 

The  plan  of  the  lecture  is  to  discuss  at  first 
the  derivation  of  the  TSD  equations  (Sec  2)  as 
it  appears  for  slender  objects  and  for  three  dimen¬ 
sional  wings.  Then  various  examples  are  studied 
where  subsidiary  expansions  are  necessary  either 
at  the  outset  or  a  later  stage.  These  examples 
are  Slender  Configurations  (Sec  3),  Lifting  Line 
theory  for  Subsonic  Flow  (Sec  k) ,  Wind  Tunnel 
Corrections  (Sec  5),  Perturbations  of  Sonic  Flow 
(Sec  6),  and  Perturbations  of  Shock  Free  Flow, 

(Sec  7).  HaturaUy  a  complete  discussion  of  these 
problems  is  not  feasible.  However  I  hope  that  the 
details  will  he  sufficient  to  allow  the  audience 
to  appreciate  how  the  judicious  use  of  asymptotic 
analysis  leads  to  an  understanding  of  the  struc¬ 
ture  of  problems,  how  it  simplifies  the  numerical 
calculations  that  must  he  done  and  extends  their 
usefulness . 


2.  Transonic  Small  Disturbance  Equations 


The  systematic  derivation  of  transonic  small 
disturbance  equations  it  discussed  in  (2],  [3],  [U], 
including  some  extensions  to  higher  order  approxima¬ 
tions.  Here  the  starting  point  Is  the  full  poten¬ 
tial  equation  which  is  valid  to  second-order  in  a 
small  parameter  for  Mach  numbers  close  to  one,  in¬ 
cluding  the  effects  of  shock  waves.  That  is,  the 
vortldty  Introduced  by  shock  waves  is  negligible. 
Thus,  within  the  gas-dynamic  framework  (inviscid, 
perfect  gas)  which  has  proved  so  useful  for  aero¬ 
dynamics  problems  the  exact  equation  for  the  poten¬ 
tial  is 


V2  *  -  9  •  v  (V)  (2.L) 


where  a  =  ■J-yR'T  -  local  speed  of  sound, 
q  «  local  velocity  «  to 


(2.1)  is  a  version  of  the  continuity  equation 
(V*  pq  =  0).  The  conservation  of  total  enthalpy  is 


a*2  7  +  1 

2  7  -  1 


(2.2) 


U  «  flow  speed  at  upstream  infinity  along  x-axls. 
a*  *  critical  speed  ■  flow  speed  and  sound  speed  at 
local  Mach  nvsnber  one.  The  local  sound  speed  is 
also  related  to  the  local  pressure  and  density  by 
the  isentropic  formulas 


a2 

2 

a_ 


p 


(2-3) 


valid  to  second  order. 


Perturbations  in  a  uniform  nearly  sonic  flow 
are  characterized  by  the  facts  that  (i)  the  stream- 
wise  perturbation  is  an  order  of  magnitude  larger 
(formally)  then  the  transverse  perturbation  and 
(li)  the  disturbances  die  out  more  slowly  in  a 
transverse  direction  then  in  a  streamwise  direction. 
Specifically  a  suitable  limit  process  for  describ¬ 
ing  the  principal  regions  of  transonic  small-di stun- 
bance  flow  has 


dimensionless 

coordinates 


and 


transonic 

similarity 

parameter 


{x  -  |  ,  y 
f  1-M2 

{-  Vs 


*=  £  I  ’  2  *  E  I  ^ 


fixed  as  6 


T 

I 


0  where  e(i),  u(: )  -0 


L  m  characteristic  length  in  the  streamwise  direc¬ 
tion,  T  »  characteristic  transverse  length  so  that 
&  «  measure  of  the  flow  deflection. 


The  form  of  the  expansion  is  free-stream  plus 
a  small  perturbation 
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♦  «  0{X  +  l  M  *(x,y,t;g)  +  •••}  (2-5) 


The  non-linear  transonic  equation  results  la  the 
distinguished  esse  that  In  »  e^l 

ED  eqn  :  (K-  (y  ♦  ♦  S2*  ■  0  , 


(2-6) 


The  transonic  coordinate  scaling  tjtua  nreser- 
ves  local  Mach  wave  structure  J*»  ey  •vny-VMjJ  -  l'y 
Tor  K  fixed.  The  local  Mach  masher  H.  is  given 
by 


«f-l 

“-■(7  +  l)*x-K  (2.7) 


so  that  the  ED  equation  being  essentially  non¬ 
linear  has  a  chance  to  represent  flow  which  is 
locally  subsonic  or  supersonic .  It  can  also  be 
expressed  in  a  conservation  fora,  (the  continuity 
equation  “  °) 

<2-e> 

whose  integrated  form  holds  across  shock  Jixnps . 
The  mass  flux  vector  is  in  this  approximation 


In  addition  vortex  drag  can  also  appear  [2] . 

m 

Dy  -  P„  uVu2(6)  JJ  \ 

(2.15) 

Finally  the  relationship  between  u  and  6 
is  found  fran  the  boundary  condition  of  tangent 
flow  at  the  surface.  Sonetiaes  these  relationships 
are  clearest  in  terns  of  an  inner  expansion  valid 
near  the  body  surface,  as  in  the  next  sections. 
Asymptotic  watching  is  used  to  relate  the  outer 
expansion  (2.5)  and  the  inner  one. 

Further  boundary  conditions  needed  to  specify 
the  flow  uniquely  are  vjpishing  of  perturbations  at 
upstream  infinity  (♦  ,V  ♦  -  0)  and  a  Kutta  condi¬ 
tion  at  sharp  trailing  edges  for  small  perturba¬ 
tions  this  latter  condition  translates  into  no  per¬ 
turbation  pressure  load  at  the  trailing  edge 
surface . 


3.  Transonic  Slender  Body. 

A  slender  body  is  characterised  by  the  fact 
that  in  (x,y,f)  coordinates  the  body  surface 
shrinks  to  a  segment  of  the  x-axis  (0  <  x  <  1)  as 
the  outer  limit  (6  -  0,  x,y,z,K  fixed)  is  con¬ 
sidered.  Such  a  body  slight  be  represented  by 


Note  that  the  streamwise  mass-flux  has  a  maximum 
where  (y  +  1)*  «  K,  or  the  local  speed  is 

exactly  sonic. 

This  is  the  approximate  three-dimensional 
representation  of  the  sonic  throat  in  a  one-dimen¬ 
sional  stream  tii>e.  The  corresponding  shock  Jump 
conditions  Included  in  (2.8)  can  be  abbreviated  as 
follfws.  let  g(x,y,z)  >0  be  the  shock  surface 
and  n  the  normal  to  it  say, 

2  »  +  Then  from  (2.9)  or  (2.6) 

across  a  shock  we  have 


B(x,r,B)  »  0  «  r-6F(x,@)  , 

0  <  x  <  1  r  «  v/y^+z2  6  -  tan-1  * 

(cf.  Fig.  3-1) 


(3-D 


In  physical  units  the  body  has  length  L.  The 
effects  of  angle  of  attack  a,  assumed  to  be  0(b), 
are  considered  here  to  be  included  in  F(x,0)  but 
can  be  written  out  more  explicitly  if  necessary.  A 
similarity  parameter  A  •  a/6  then  appears . 


The  boundary  condition  of  tangent  flow  on  the 
surface  q  ■  VB  >0  is  thus 


«x[»x-  H*  *x]  +  E^ls  •  0 

Also  from  continuity  of  tangential  velocity 

[•],-0  (2.11) 

where  the  Jump  of  a  quantity  at  a  shock  [  )  - 

(  ),  -  (  )-,  (  ).  -  quantity  evaluated  behindthe 
shock  (  )°  «  quantity  evaluated  ahead  of  the 

shock. 

Consistent  with  the  approximations  made  the 
usual  pressure  coefficient  is  given  by 

cp  -  -2u(6)*x  (2.12) 

Further  an  expression  for  the  wave  drag  Du 
can  be  derived  within  the  XSD  framework  [2]  * 

-  -P„  AV(»)  ^ JJ  t^]3  dydl  (2.13) 

shocks 

This  has  a  simple  physical  interpretation  as  the 
integrated  entropy  Jmp  across  the  shock  waves 
since  for  transonic  shock  waves 


*r  "  wx  Fx  +  "I  ®e  Fe  00  r  -  BF(X>6)  (3-2) 

Because  of  the  singularities  in  $  as  r  -  0  an 
inner  coordinate  r*  »  r/6  is  useful  in  which  the 
scale  of  the  body  is  preserved  r*  ■  F(x,0)  on  the 
body  surface.  The  corresponding  limit  process  has 
6-0  (x,r*,6;K)  fixed.  The  inner  expansion  is 

of  the  form 

•(x,r,0)  •U{X  +  L(T0(6)({!0(x,r*,e)  (3.3) 

♦  tj(6)<Pl(x,r*,e)  +  T2(6)*2(x,r*,e))+-  •  •) 

t>0  turns  out  to  be  a  switchback  term  introduced 
for  purposes  of  matching  and  to  depend  only  on  x 
(see  below). twit chback  is  due  to  the  occurrence  of 
log  terms  in  which  is  the  term  connected  with 
the  boundary  cOiditloos  (3*2).  We  have  from  the 
equation  for  the  exact  potential  (2.1) 


*  * 

°  ar*2  7* 

since 


*0  1  A0 

V*2  <F1  -  0 


(3-1*) 


,3  lii 


ET 


2>(a  ]3. 

1  x's 


The  inner  flow  is  a  solution  of  a  cross-plane 
laplacian .  The  flow  is  analogous  to  unsteady  in¬ 
compressible  flat  in  (r*,9)  with  x  appearing  as 


2 


tine.  The  boundary  condition  (3.2)  expressed  in 

inner  variables  is 

t,(£) 

— j U.F.e)  ♦  •••  ■  (3.! 


6{1  4  •••)Fx  4 


^(6)6  d>le(x,r,e)p6 


This  shows  that  T^b)  ■  6Z  and 

*1  #  “  Tx  *  V’lg  ®  r#  "  (3 

The  order  of  magnitude  of  the  perturbation  poten¬ 
tial  near  the  body  is  thus  proportional  to  the 
ratio  of  cross-section  area  to  square  of  body 
length.  Note  that,  on  the  surface, 

»1 

jn*  *  D  -^1° 

n*  *  normal  to  surface  in  a  cross-plane  ■  'i  , 
B*  -  r*-F(x,e)  .  0 


i  -I  —  F 
*r*  *e  r* 


on  r*  «  F(x,e). 


S)} 


*  pa2  e 
on  r*  •  F(xpp) 


1  *  ^*2  Fe 


from  (3-6).  Now  asymptotic  matching  of  the  inner 
and  outer  expansions  in  this  case  is  concerned 
with  the  behaviour  of  the  inner  expansion  as 
r*  -  •  and  the  outer  expansion  as  r  -•  0 . 

A  class  of  intermediate  limits  has 

(x,r  ,e,K)  fixed 

T1 

as  t  -  0  where  ^ 

r,i  ”  ^7  *uch  that  K2(fc)  » ’i  >>  6 


Thus  r*  »  - 


•  2  2 
r  ■  u  r  ■  u  r,r  -  0  . 


The  corresponding  asymptotic  behaviour  of  inner  and 
outer  solutions  must  be  calculated  to  carry  out  the 
matching.  For  the  inner  solutions  <?.  we  have  the 
integral  theorem  L 

n-  /(•£»  -’is?)"*  (3-9) 

where  1 

0  “4  log((y*-y*)2  4  (,*-«*)2}2 

(y|,i£)  »  coordinates  in  surface. 

which  expresses  the  potential  as  a  distribution  of 
sources  and  doublets  cn  the  surface  (cf.  Fig  3.2) 
Thus  as  r*  -  •  from  (3.9) 


*1  -  aF log  r# 

or  from  (3.7) 


H»): 


/r2  4  if  +  •■•(3.1fl) 


-  Sr  *•«  '*/' 


FFxd6  4 


-  ^T^log  r»  4  g*<x) 
where  A(x)  -  \  J^F^x.e) de 


(3-11) 


-  cross-section  area  (dimensionless) 

The  effective  cross-plane  source  strength  is 
given  by  the  streasiwiae  rate  of  change  of  cross- 
section  area.  Further  terms  in  (3.U)  can  include 
g*(x)  and  terms  which  die  out  as  r*  -  ••  The 
outer  equation  (2.6)  In  cylindrical  coordinates  is 


(X(^>\>Vxx 


♦♦rr+*  Vp>*ee*°  (3-12) 


but  for  matching  no  singularity  of  larger  order 
than  log  r  should  appear  as  r  -•  0.  The  solu¬ 
tion  of  (3.12)  must  behave  near  the  axis  as 


•(x,r)  -  B(x)  log  r  4  g(x;K)  4 


0<x<l  (3.13) 


Now  to  carry  out  the  matching  of  inner  and  outer 
expansions  under  the  Intermediate  limit  we  must 
have 

outer 


u(&KE(x)  lo^i2T,rJ  4  g(x)  4  •••l  4-* 

^  Inner  ^ 

t0(6)Vo(*)  +  lo«  +  + 


It  follows  that 


■<*>  * 


(3.13) 


g(x;K)  »  g*(*»K)  (3-1*0 


and  the  switch  back  tens  can  be  chosen  Tn(t)  = 
&  log  52  ,,  0 


(3-15) 


The  situation  of  the  dominant  term  is  new 
clear.  Whether  or  not  the  slender  body  is  axi- sym¬ 
metric  or  lifting  the  dominant  outer  flow  *(x,r) 
is  axi- symmetric.  The  first  term  of  (3.I3)  serves 
as  a  known  boundary  condition  which  allows  *(x,r) 
to  be  calculated  as  a  solution  of  the  non-linear 
transonic  equation  following  the  numerical  methods 
for  example  of  [3].  This  determines  the  unknown 
g(x).  For  further  details  see  [6]. 

Since  the  outer  flow  is  completely  specified 
by  the  cross-section  area  distribution  of  the  body 
any  shock  waves  which  occur  are  also  related  only 
to  this  quantity.  The  wave  drag  Eu  of  (2.13)  is 
expressed  as  the  Integral  of  the  entropy  jump  over 
the  shocks  and  thus  depends  only  on  the  cross- 
section  area  distribution.  The  wave  drag  is  that 
of  an  equivalent  body  of  revolution  which  has  the 
same  axial  distribution  of  cross-section  area 
(Transonic  Equivalence  Rule).  Chan  [7]  has  inves¬ 
tigated  a  special  class  of  bodies  of  revolution  to 
see  which  has  the  minimus  wave  drag  for  a  given 


3 


(K,E).  However  it  eeeiu  likely  that  ebook-free 
axially  symmetric  chapes  axlet  which  bare  aero 
wave  dr as  for  high  subsonic  epeeda,  although  none 
have  been  calcinated  yet. 

The  procedure  outlined  here  can  be  carried  to 
higher  order  to  obtain  correction*  to  the  domi¬ 
nant  terse  considered  here.  In  eussary  the  foes 
of  the  Blender  body  expansion  for  the  potential  la 

•  -  UL{x  ♦  62*(x,?)  +  0(6^)}  outer 

-  UL  jx  ♦  62  log 

♦  62$j(x,r*,e)  +  0(6**  log2  6  )|  inner 

The  switchback  tars  for  tranaonic  flow  ia  larger 
than  that  for  Unearned  aubaonlc  or  auperacnic 
flow. 


The  treatsent  of  lifting  line  theory  ia  a 
little  different.  The  atarting  point  ia  a  IBS 
formulation  of  lifting  aurface  theory.  The  system 
ia  planar  in  the  sense  that  aa  6-0  the  body 
aurface  collapses  to  a  plane  ^  -  0.  For  this 
ease  we  can  expect  ISO  theory  to  be  valid  froa 
the  aurface  to  infinity  with  the  following  under¬ 
standing.  There  la  a  nan -uniformity  near  infinity 
with  the  unrolled  trailing  vortex  sheet  located  at 
y  ■  0,  and  local  nan-unifoxsitiea  near  stagnation 
points.  These  non-unifoaraities  do  not  affect  a 
good  first  approximation  to  surfnae  pressure  dis¬ 
tribution,  lift  and  drag.  The  dependence  of  the 
TSD  solution  on  scaled  aspect  ratio  is  studied 
for  large  aspect  ratio.  The  first  approxlmtlon 
near  the  wing  Is  the  flow  for  Infinite  aspect 
ratio  (two-dimensional  f  Imf  and  tbs  second  term 
represents  a  correction  to  this  flow.  The  flow 
far  fros  the  wing  baa  as  an  epproxlamtion  the 
everywhere  subsonic  flow  past  a  vortex  line  with 
trailing  vortex  sheet .  The  near  field  and  far 
field  flow*  natch  asymptotically .  This  approach 
to  lifting  line  theory,  which  unlike  Prendtl's 
does  not  require  the  solution  of  an  lntagral  aqua¬ 
tion,  was  pioneered  by  VanDyks  for  the  incompres¬ 
sible  case. (8].  The  detailed  application  of  these 
ideas  to  transonic  flow  ia  given  in  [9]  • 


The  characteristic  length  L  is  chosen  as  the 
airfoil  chord  and  the  upper  and  lower  surface  are 
represented  by  Fu  gi  b  Is  the  (dinensionless) 
span,  (cf.  Fig.  k.i; 

W(x,y,«)  -0  -y-wu>J(«,  £)♦«*  Saw*. 

r  (4.1) 

The  representation  here  regards  Ftt  /  as  a 
fixed  surface  and  Q  is  the  angle  of  gitack.  The 
boundary  condition  of  tangent  flow  (q  •  W  ■  0) 
can  be  applied  approximately  at  J  ■  0  and 
becomes 

u3/2  ♦yfx.O*)*  •••-6  -a,  B-b61/3.  (4.2) 

Thus  for  a  distinguished  ease  u3i  2  ■  6  and 
6  -  O(o)  or  Or/6  ■  A  fixed  in  limit  as  6-0 
is  a  similarity  parameter.  The  TBS  mansion  is 
thus  associated  with  the  limit  6-0  (x,y,f  ;K,A,B) 
fixed  and  has  the  form 


•  -  UL{x  +  62/3*(x,M;K,A,B)  +  6U/3*(2)+  •••  (4.3) 

(of.  Fig.  4.2).  The  leading  and  trailing  edges  are 
x  *  x^  t(2/B).  The  Kutta  condition,  by  (2.12),  is 

■  Vv0+,i)  ‘Vv0-*^  “  0  - 

The  lift  Z  is  given  by  the  spanwise  integral  of 
the  dreulatian  distribution  T(t). 

t  -  Tii)dz  (4.5s) 

where  -B 


r<*> 


(4.5b) 


1 J**  (*x(x,0+,I)  -  *x(x,0-,t))dx«i  f 


Across  the  trailing  vortex  sheet  there  is  no  Jump 
in  pressure  so  that 

t**]^  ■  *x(x,0+,l)  -  *x(x,0-,5)  «0  x  >xT  (4 .6) 


Bow  the  dependence  of  ♦  on  aspect  ratio  B  is 
considered  for  B  -  •.  In  the  outer  limit 

(x*  -  y*  -  |,s*-|jare  fixed  as  B  -  »  ; 

the  plsnform  shrinks  to  s  line  (of  singularities). 
Transverse  length  scales  are  measured  essentially 
in  terms  Of  .the  span.  In  the  inner  limit 


Stems  of  the  span. 

,7,t*  ■  1)  are  fixed  as  B  -  the  relative 
amrise  location  is  held  fixed.  The  resulting 
forms  of  expansion  are 

Inner:  (4.8) 

*  -  ♦()(*  »y»**)  +  |  •1(x,?;z*)  +  ••• 

Outer:  (4.9) 

♦  -  <J>0(x*,y*,s*)  <Px(x*,y*,**)  * 

+  5  5P2(x*,y*,x*)  +  ••• 

♦0  represents  the  two-dimensional  transonic  flow 
past  the  airfoil,  including  shock  waves  if  they 
occur.  The  aspect  ratio  correction  is  0(1/B). 

♦.  satisfies  the  vmriatlonal  equation  of  two- 
dimensional  transonic  flew 


(K-(y  +  !)♦.)♦,  -  (y  +  1)  +  ♦, 


.0  (4.10) 


This  equation  is  linear  but  of  mixed  elliptic 
hyperbolic  type.  The  variable  coefficients 
♦Q  ,  *0  are  regarded  as  known  and  the  elliptic 
x  uxx 

and  hyperbolic  regions  for  a  given  solutions  are 
the  same  as  that  for  «...  Since  *.  takes  account  of 
the  airfoil  shape  the  correction1  airfoil  looks 
like  a  flat  plats 

♦  ^(x, (>♦;**)  »  0  xJ((**)<x<xT(z*)  (4.11) 

Iftx*T*  “  0  and  ( *^)rg  *  represents  the 
effect  of  the  Kutta  condition. 

In  addition  there  are  shock  conditions  which 
are  calculated  by  perturbing  the  shock  location. 


k 


*  -  e^i**)  ♦  5  g^y.**)  ♦  •••  (^-12) 

and  expanding  about  the  original  chock  location. 

In  this  first  treatment  the  very  weak  singularity 
at  the  foot  of  the  shock  is  neglected,  a  better 
theory  should  use  local  equations  near  the  shock 
analogous  to  the  development  in  Sec.  7  below 
(cf.  also  [13]). 

Numerical  calculations  can  be  carried  out  to 
find  the  nonlinfear  two-dimensional  potential  *0 
including  shock  waves  by  one  of  the  standard 
methods.  This  provides  the  first  approximation 
solution  as  well  as  the  variable  coefficients  for 
the  variational  equation  (4.10).  The  far  field 
boundary  condition  for  the  variational  equation 
comes  from  matching. 

In  the  outer  expansions  the  term  is 

a  switchback  term  introduced  for  purposes  of  match¬ 
ing.  The  equations  valid  at  some  distance  from  the 
wing  are 

+  Vn  *  4n  -  0  (4.X3) 


0  (4.14) 

(>  +  l)<P0  (4.15) 

x*  x*x* 


is  given  by  linearised  flow  past  a  vortex-line 
with  trailing  sheet 

s-5  r-r^f  =£=W 

4'  J-1  y**>(z*-i)  l 


♦ 

40 

+  4, 

°x*x* 

y*y 

1 

*1 

■f 

*  <p. 

y*y* 

«P2 

•f 

4> 

+  4. 

x*x* 

2y*y* 

i 

(4.16) 

The  corresponding  representations  for  cp^  g 

involve  higher  singularities  [9] .  Matching  is 
carried  by  a  class  of  intermediate  limits  where 


re  * 


is  fixed  es  B  -  »  ,  P  -  » 


?  ■  vx2  +  $2  »B r*.  Thus  r  «  £( B)r ^  -•  •  , 
r*  * !  rP  -  ° 

Then  the  far  field  of  the  inner  solution  in  each 
cross-section  plane  of  the  wing  matches  to  the 
singular  behavior  of  the  outer  solution  as  the 
lifting  line  is  approached.  For  example 

2 

,  r0n  .  ly  +  UTo\loe.  T  .  „  .  if  D0  _  „ 

*0  *  -  F  9  *  I  K\^j“t"C02Se+?t^K 

+  rr: 6in  6  -  i  )COi  + '  •  ’  > ? — 

2rVK  (4.17) 

V—T=“e’SJ1  tttt*!  ♦  *  (**.y*)— 

-1  (4.16) 

are  used.  Matching  shows  that  the  circulation 
around  the  two-dimensional  wing  is  that  around  the 
vortex-line  yQ  »  rQ  and  further  that  the  varia¬ 
tional  equation  potential  has  an  induced  downwash 
at  infinity  as  a  result  of  the  vortex  sheet. 

5  r1  r6<o  rie 

*l>x#,y>x*)  m“  t-rj  j  irrjdi-TT+”  M  ■ 

(4.19) 


(4.19)  is  the  missing  boundary  condition 
necessary  to  complete  the  problem  for  in  (4.10). 

R.D.  Small  [10) ,  [11]  has  devised  the  necee- 
sary  type-sensitive  algorithm  for  (4.10)  and  a 
special  way  of  treating  the  shock  which  appears 
For  the  case  of  similar  airfoil  sections  calcula¬ 
tions  need  to  he  carried  out  in  one  cross-section 
plane  only.  Sac  of  Small's  results  for  a  elliptic 
plan  form  with  NACA  0012  sections  appear  in  Fig. 
4.3.  The  shift  of  the  shock  location  is  clearly 
•een. 

The  approach  initiated  here  has  been  carried 
further  by  Cook  [12]  to  treat  the  case  of  swept- 
back  wings  where  additional  log  terms  are  needed  in 
the  inner  ex pans i on .  Numerical  calculations  and 
some  further  analysis  based  on  Cook's  ideas  appear 
in  the  paper  of  Cheng  et  al.  [13] . 

The  formulation  of  lifting  surface  theory  at 
the  beginning  of  this  section  can  also  be  used  as  a 
starting  point  for  slender  wing  theory.  Pointed 
planforms  can  he  considered  under  the  limit  B-0. 
The  inner  expansion  leads  to  Laplace's  equation 
in  cross-section  planes  but  the  primary  outer 
expansions  has  a  dominant  term  which  satisfies  the 
non-linear  transonic  equation  (2.6).  The  dominant 
reaulte  of  the  previous  section  are  retrieved  for 
the  case  of  slender  planar  planforms.  The  expan¬ 
sions  have  the  form: 


♦  -  B  log  B3/2<^(x)  +  Bq0(x,y*,z*)  + 

+  BU  log  B3^2q.  + 

+  B1*  log  B3/2<f2  +bVj+  •  •  • 

(4.20) 

Outer 

♦  «  B  *0(x,y,z)  +B3/2*1(x,y,5)  +  •••  (4.21) 

where  y  «  -/b  y,  z.  *  Vb  z.  It  also  is  necessary 
for  these  limits  to  make  sure  that  the  original 
K  -  0  such  that  1  -  M2 

R  -  |  .  bg  "  is  fixed  .  (4.22) 

Corrections  to  the  slender  wing  lift  and  drag  can 
he  found  this  way. 

5.  Subsonic  Wind  Tunnel  Effects. 

Wind  tunnel  correction  effects  can  be  calcu¬ 
lated  by  following  a  plan  similar  to  that  of  the 
previous  section.  If  the  tunnel  height  is  con¬ 
sidered  large  the  first  (inner)  term  is  transonic 
flow  past  the  wing  in  a  free-field  and  the  second 
term  represents  a  correction  due  to  the  presence  of 
walls.  This  correction  is  evidently  connected  with 
a  solution  of  the  variational  equation.  Near  the 
tunnel  walls  however  the  flew  has  a  structure 
highly  dependent  on  the  wall  boundary  conditions. 

A  separate  outer  expansion  is  needed  to  describe 
this  region.  Matching  of  the  inner  and  outer 
equations  provides  the  necessary  far-field  boundary 
condition  for  the  solution  of  the  variational 
equation. 

The  TSD  problem  for  a  two-dimensional  air¬ 
foil  (L  ■  chord)  in  wind  tunnel  with  solid  walls 
is  considered.  Following  the  previous  section  the 
ISC  expansion  has  the  form 
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*  -  U^fx+B2/3  ♦(x.y}H,l^,AT)  +  •••  }  (5-1) 

where  •  upatream  velocity  la  the  tunnel 

Hj  •  upstream  Mach  number  in  the  tunnel 

o_ 

•  -gS  Oj  -  angle  of  attack  in  the  tunnel 

H  ■  h6^3  ■  tunnel  aimilarity  parameter, 

!  tunnel  height 
“  “  chord 

(cf.  Fig.  5-1) 

The  inner  expansion  la  connected  with  the  lialt 
H  -  ■  (x,y)  fixed.  Further  since  corrections  are 
being  sought  a  relationship  between  tunnel  and 
free-flight  parameters  is  assumed 


Hatching  is  carried  out  with  a  intermediate  limit 
r 


fixed  0«n«P 


such  that 


r.  .  4S1 

as  r*  -•  0  we  note  from  (5. 10) 

9, 


r*  .  ^  ~  0  (5 .11) 


12) 


while  hsi  other  singular  terms  >  J"; 

not  reproduced  in  detail  here.  The  far  field 
of  the  inner  potential  correction  has  the  form 

r 

*1  -  BixS'4Elx+M1y*2f  ®  +  (5’13) 


^  -  X,  ♦  U^H)^,  At-A^  +  u1(H)Ac  . 

(5-2) 

The  inner  expansion  valid  near  the  airfoil  is  of 
the  form 

♦  -  ♦0(x,5;Kp,Ay)  +  u1(H)*1  +  •••  (5-3) 


which  leads  to  the  TSD  equation  and  a  variational 
equation  modified  for  correction  parameters 


< V  +  *1~  “ 0 

(L-Oy+D*  )♦.  -  (y+l)*.  ♦ 

X  XX  X 


-  0 


(5-4) 

(5-5) 


Matching  shows 


*0 

Mi-5? 


ui(»)  -  4 


r_  ? 

B  .  -I  -L 
1  ^  2k 


V 


2  h-Jkf 


where 


]wto'rF 


The  x5r  in  (5-13)  matches  to  the  x«y*  in 
(5*12)  which  is  the  near  field  of  a  point  vortex 
reflected  in  the  walls.  The  upwash  term  M.y 
matches  to  the  reflection  of  a  divortex  in  walls 
part  of  9  . 


The  outer  expansion  preserves  the  tunnel  geometry 
so  that  the  airfoil  shrinks  to  a  point  (singular)  as 

H  -  •».  The  limit  process  has  x*«jj  *  y**{j  fixed. 
The  outer  expansion  is  of  the  form 


♦  -  <P0(x*,y*)  +  *°§  H  ^(x^y*)  + 


where  the  switchback  term  0 
introduced  for  matching.  The' 


Vo  .  *  ♦  *0  .  * 

J  ^  v“v“ 


Vv. 


.  *  *1 


y#y* 


(7+l)<Pn  % 


0  ”  Vo 

x*x*  x*x* 


(5.9) 


(cf.  Fig.  5-2). 

The  far  field  for  the  first  inner  potential  is 
(1.17)  with  X  -  K_,r  r_,  D.  -  D0_  etc.  Thus  for 
r0  the  airfoil  appears  as  a  vortex' and  the  non¬ 
linear  terms  in  (5*9)  farce  e  modified  behavior  of 
higher  order  terms  which  provide  a  far-field  for 
He  have 

*o  +  i*o  *  5T  lo«  * 


Z*  •  X*  ♦  i 


/Xyy* 


r*e 


16 


(5.10) 


In  using  this  method  a  free-field  calculation 
for  *0  is  made  and  then  for  some  choice  of  cor¬ 
rections  (K-tA  )  the  variational  equation  is 
solved  and  a  wind-tunnel  correction  field  is  found. 
It  is  not  possible  in  general  to  find  a  combination 
of  parameters  (XC,A  )  such  that  the  wind  tunnel 
results  reproduce  free  flight  pressure  distribu¬ 
tions;  some  minimization  could  be  performed - 
However  it  is  possible  to  choose  (K  ,A  )  such 
’that  the  C.  «  C.  .  That  is  angle  8f  attack  and 

rr  r 

X.' .uaber  corrections  are  fbund. 

he  tails  of  the  method  outlined  here  appear  in 
[14].  The  method  was  applied  earlier  by  Chan  [15] 
for  porous  walls  but  there  is  seme  difficulty  in 
reconciling  his  results  to  the  case  considered 
here .  A  curve  showing  lift  corrections  against 
A  (X  »  0)  which  is  a  result  of  preliminary 
numerical  calculations  is  shown  in  Fig.  5.3. 

6.  Perturbation  of  Sonic  Flow. 

Perturbations  of  exactly  sonic  flows  for  one 
reason  or  another  also  involve  non-uniformities . 

The  case  where  M  -  1  from  below  is  particularly 
clear.  For  M  <  1  the  transonic  far-field  of  a 
two-dimensional  airfoil  is  basically  subsonic  and 
given  by  the  form  (4.17).  Frcm_this  form  it  can 
be  noted  that  as  X  -  0  (M  -  l")  larger  and 
larger  distance!  7  from  fhe  origin  must  be  con¬ 
sidered  if  this  description  is  to  be  valid-  For 
1*0  (M  •  1)  this  far  field  disappears  and  is 
replaced  by  a  similarity  form  characteristic  of 
X  «  0,  which  incidentally  always  has  a  potential 
♦  ayntric  about  5*0.  Thus  if  we  try  to 
represent  the  perturbation  flow  *  for  small  K 
as  that  at  aocic  plus  a  small  correction  there  is 


Nr 


a  non -uniformity  at  infinity.  However  near  the 
body  this  representation  is  valid  and  leads  to  a 
"lav  of  stabilization"  or  Mach  nvmber  "freeze". 

It  is  of  interest  to  obtain  the  order  off  the  cor¬ 
rections  and  a  way  to  calculate  them.  Another 
problem  in  which  two-dimensional  sonic  flow  is 
perturbed  is  the  finite  aspect  ratio  wing  at 
M  >1.  The  spirit  of  the  approach  used  for  this 
problem  is  that  of  the  lifting  line  theory.  A 
brief  discussion  is  given  of  both  of  these 
problems . 

The  nearly  sonic  flow  past  a  two -diaens 1  cue 1 
airfoil  was  recently  analyzed  from  the  point  of 
view  presented  here  by  Cook  [16].  Some  earlier 
Russian  papers  [17],  [18],  [19]  treat  the  same 
problem  without  a  clear  discussion  off  asymptotic 
matching  and  validity. 

Central  to  the  discussion  is  the  similarity 
solution  representing  the  asymptotic  far  field 
past  an  airfoil  at  -  1  and  its  possible  per¬ 
turbations.  The  general  structure  of  this  flow 
is  shown  in  Fig.  6.1.  The  body  appears  as  a 
point  (singular).  The  far  field  has  sonic  line, 
limiting  Mach  line  (characteristic)  and  shock  wave 
all  lying  cn  similarity  curves.  The  flow  up  to 
the  limit  characteristic  is  found  as  a  mixed 
elliptic-hyperbolic  problem  and  then  the  flow 
after  the  limit  characteristic  and  the  shock  is 
found  by  analytic  continuation.  The  flow  down¬ 
stream  of  the  shock  is  supersonic.  The  asymptotic 
form  is  2  ^ 

•(*,y) -j^a^ffa?) +cQ  +  ciy  5  a"3fx(ai)  + 

(6.1) 

vber-- 

5  “  T* 

f  satisfies  a  non-linear  ordinary  differential 
equation  and  f^  a  corresponding  linear  varia¬ 
tional  one.  The  closed  form  solutions  for  f,f. 
have  been  worked  out  by  various  authors  of  which 
we  cite  here  only  [20].  The  principal  boundary 
conditions  are  regularity  on  J  ■  0  and  smooth 
passage  through  the  limit  Mach  wave  ■ 

The  particular  exponent  2/5  and  similarity 
curve  {  »  =  4/5  is  isolated  ttm  the  general 

„  J 

form  jr  f(|),  |  •-  by  these  conditions. 

f 

Further,  smooth  homogeneous  similarity  solutions 
of  the  sonic  variational  equation  of  the  form 

yVe)  exist  only  for  o  «  |  n,  |  ♦  |  a,  -|n  -  j  , 

r n  ♦  1  where  n  >  integer.  These  considerations 
ire  important  for  matching  and  are  discussed  in 
[20],  [22]. 


f..3A  .-5(^:|.,6),  ,..5  .-5(2^1)  , 

-  AV5  <6's) 

where 

0  <  a  <  4/3  “  limit  Mach  line  . 

The  limit  process  for  the  inner  expension  in 
this  case  has  (x,J)  fixed  as  K  -  0,  so  that  the 
sonic  region  grows.  The  corresponding  form  turns 
out  to  be,  after  matching, 

♦(x,f;K)  -*0(x,y)  +K^I  +  K3/2*2(x,y)  +  •••(6.3) 

with  (6.3)  valid  near  the  airfoil. 

The  limit  process  for  the  outer  expansion  must 
grow  as  fast  as  the  sonic  zone  and  preserve  the  non¬ 
linear  structure  of  the  dominant  equation.  Also  the 

form  of  similarity  curves  ^  should  be  pre¬ 

served.  The  limit  process  has  (x,y)  fixed  where 

x  -  K2  x  ,  y  »  (6.4) 

and  the  corresponding  outer  expansion  is 

♦U.y.iU  -  t  ?0(x>y)  ♦  <M*,y) ♦  K^2M*,y)  +  •  • 

K  0  1  2  (6.5) 

The  boundary  value  problems  to  be  solved  for 
the  inner  appr  cad  mat  lens  are 

-(Y*l)  ♦„  ♦„  +  *0  -  0  (6-6) 

x  xx  jFy 

♦«  (*>&+)  ■  C  #(*)  “  the  airfoil  (6.7) 

Uy  —  U,* 

♦0U,y)  -y^a'3  f (ai)  +  c0  +  c^y  V3fx(a{)  + 


should  be  pre- 


Tbat  is,  ♦q  is  the  potential  of  exactly  sonic 
flat  past  the  airfoil. 

A  class  of  flows  past  realistic  airfoil  shapes 
at  Hm  m  1  has  been  calculated  by  S.  Tse  [23] . 

Tse  used  a  hodograph  representation  and  solved 
Triccmi's  equation  nimerically  using  type-sensitive 
differencing  with  the  proper  free-stream  singularity 
corresponding  to  (6.8).  Tse's  solutions  for  planar 
flow  are  also  used  by  Ziegler  [22]  for  sections  of 
a  three-dimensional  wing.  These  solutions  can  be 
considered  to  give  ♦„  for  a  special  class  of 
shapes.  A  typical  holograph  and  airfoil  shape 
are  given  in  Fig.  (6.2),  (6.3),  (6.4). 

The  principal  correction  potential  satisfies 
the  variational  equation 

-(y+l)(c.  ♦.  )  ♦  *  -  0  (6-9) 

xx 


Germain  [21],  by  constructing  an  integral 
theorem,  showed  how  to  relate  the  scale  factor  _s 
of  the  far  field  to  the  solutions  on  the  front  " 
part  of  the  airfoil  up  to  the  limiting  Mach  lines. 
Ziegler  [22]  has  extended  Germain's  ideas  to 
obtain  an  alternate  derivation  of  the  formula  for 
a  and  also  for  c,.  One  parametric  representa¬ 
tion  of  the  closed  roro  is 


(x,0+)  »  0  cn  the  airfoil 

2y  ' 

and  fran  matching 

♦2  -  d2  5°  y  -  • 


(6.10) 
(6. 11) 


For  the  terms  of  the  outer  expansion  it  is 
found  that 

(l-(-yH)*0_)<f0__  +  <P0__  -  0  (t .  12) 
x  xx  yy 


7 


-  y  8/,5»‘3  f(»t)  +rrr+Dny  ♦  aa  y  - 


Outer  ttnaanslon 


*.y  -  • 


(6.13) 

(6.14) 


with  corresponding  equations  for  the  higher  term*. 
Matching  is  carried  out  with  a  suitable  class  of 
intermediate  limits  which  preserve  similarity 

curves  -j ■  j.  Matching  shows  that 

o  -  1  ,  dj  »  Iq  .  (6.15) 

The  solution  of  the  boundary  value  prcblan  for 

(6.12,13,14)  which  Includes  a  shock  wave  deter* 
nines  D_.  This  then  c expiates  the  boundary  value 
problem  for  the  inner  correct  Iks  *£  (6.9,10,11). 

The  problem  of  sonic  flow  past  a  finite  aspect 
ratio  unswept  wing  is  considered  briefly  here. 
Preliminary  results  appear  in  (24).  The  problem 
was  originally  considered  by  Ouderley  using  (non¬ 
linear)  hodograph  equations  [25].  He  was  able  to 
give  the  general  farm  of  the  corrections  and  the 
results  obtained  in  [22]  agree  with  his.  The 
problem  is  carried  further,  although  not  yet  com* 
pletely  solved.  For  this  problem  we  expect  a  sonic- 
surface  and  limiting  Mach  surface  analogous  to 
those  in  two-dimensional  flow.  Hsnce  the  mixed 
flex  ahead  of  the  limiting  Mach  surface  can  be  con¬ 
sidered  first  and  then  the  supersonic  flew  behind 
the  limiting  Mach  surface  can  be  considered.  For 
the  problem  considered  the  lifting  surface  theory 
is  as  in  Sec.  4  with  K  >  0. 

Here  we  stmssarlxe  only  the  limiting  processes 
and  tie  forms  of  the  carreapaiding  expansions.  The 
first  term  of  the  inner  expansion  represents  the 
two-dimensional  sonic  flow  past  an  airfoil  section 
and  the  first  (nan-trivial)  term  of  the  outer  solu¬ 
tion  represents  flew  past  a  line  (singular)  with 
tail  shock  and  trailing  vortex  sheet.  The  vortex 
sheet,  of  course,  does  not  affect  the  flow  over  the 
wing  surface. 

The  inner  expansion  has  the  limit  process 
(x,y,s*)  fixed  as  B  —  •  where  s*  «  j  ■  jj. 

Utder  this  limit  the  relative  spanwlse  location  is 
fixed  and  planar  transonic  coordinates  are  preserved 
leading  to  the  two-dimensional  non-linear  equation 
for  the  first  tern.  In  the  outer  expansion 
(jc*>y*,x*)  are  fixed  as  B  *•  »  where 

x*  «  '  y*  *  jj*  The  order  of  the  first  term 

is  adjusted  to  guarantee  the  non-linear  TBS  equa¬ 
tion  in  these  coordinates  so  that  the  sonic  line  and 
mixed  flow  can  extend  to  infinity.  Utader  this 
limit  the  airfoil  surface  collapse  to  the  line 
x*  ■  y*  ■  0.  The  higher  terms  in  the  extension 
satisfy  forced  or  unforced  variational  equations. 

The  results  are: 


♦(x.y.x)  ■  B8/5C0(»*)  +  e0(x,y;t*)'fB  5*x(x, fit*)*' 


both  ahead  and  behind  of  shock 


(6.16) 


♦(*,y,«) -^59Q(x*,y*,f)  +<P1(x*,y*,**)  -r 

+  B*1//5q>2(x*,y*,z*)  +  ••• 


ahead  of  the  shock 


B2/5$0(x*,y*,x)  ♦  B^jfx*.}*,**) 


♦  d>2(x*,y*,**) 


(6.17) 


behind  the  shock  (6.16) 

Matching  is  carried  out  in  the  class  of  inter¬ 
mediate  limits  (c*,x  ,y.)  fixed  where 
^  P 

*,  -  ’l(B)*  ,  yp  -  P(B)y  ,  (6.19) 

and  q  >>B-1*^ 

*  y« 

y  *  -*•  . 

’i  p 

1  ‘fl5m  77§-“T^“conBt- 

The  aspect  ratio  correction  to  the  wing  pressure 
distribution  is  thus  0(B~^). 


7.  Perturbation  of  Shock  Free  Flow 

For  a  special  airfoil  shape  flying  at  a  speci¬ 
fied  dealpi  Mach  number  and  angles  of  attack  it  is 
possible  to  have  a  shock  free  siixed  flow.  Many 
such  shapes  have  been  calculated  by  hodograph 
methods  [26],  [27].  A  supersonic  region  appears 
over  the  airfoil.  It  has  been  proved  that  these 
solutions  are  isolated  with  respect  to  small 
changes  in  airfoil  shape  [26].  This  means  that 
there  are  no  smooth  neighboring  solutions  and  thus 
that  shockwaves  appear.  The  solut; on  under  design 
conditions  has  zero  drag  but  the  neighboring  solu¬ 
tions  must  have  drag.  Presixnably  the  shock- free 
solutions  are  also  Isolated  with  respect  to  changes 
in  free  stream  Mach  number  and  angle  of  attack.  It 
is  thus  of  interest  to  calculate  the  variations  in 
the  flex  around  the  design  condition  to  see  how 
rapidly  the  drag,  for  example,  changes  with  depar¬ 
ture  from  design .  Cven  more  important  is  same 
understanding  of  what  these  departures  depend  on. 

This  problem  is  not  solved  but  some  general 
remarks  can  be  made  here  about  the  formulation, 
limiting  processes ,  and  inner  and  outer  expansions . 

In  a  TBD  formulation  changes  in  the  similar¬ 
ity  parameter  from  its  design  value  K  for  shock 
free  flat  can  he  considered  a  small  parameter. 


e  -  K-Kd 

shock-free  flow  exists  for  K  ■  K . 

a 

(V^Vxx^yy -° 

♦y(x,0+)  -  F;,(x)  0<x<  1 

0  at  • 

lx 


That  is 
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with  Kutta  condition  baa  a  smooth  aolution. 

The  outer  expansion  valid  away  from  the  rear 
of  the  smooth  supersonic  sene  where  a  shock  is 
expected  to  form  is  connected  with  the  limit  pro¬ 
cess  (x,?)  fixed  as  c  -  0.  It  has  the  form 

♦(x,y;K)  -  ♦0(x,y;Kd)  +  e  ♦  •••  (7-5) 

where  *0  is  the  smooth  solution  of  (7 >2, 7 *3 >7 ^ 
•1  satisfies  the  variational  equation 

(K  -(-y+l)*  )^-(Ytl)*0  -*0  (7-6) 

xx  x  yy  xx 

This  linear  equation  can  not  allow  shock  waves  and 
so  cannot  be  valid  in  the  neighborhood  of  (x  •  xd> 
y  ■  0)  where  the  rear  part  of  the  sonic  line 
intersects  the  surface  (cf.  Tig.  7.1).  A  small 
shock  is  expected  to  develop  in  this  neighborhood. 
In  order  to  describe  the  formation  of  this  shock  a 
local  non-linear  equation  has  to  be  derived.  A 
suitable  limit  process  has  (x+,y*)  fixed  as 
e  -  0  where  _ 

**  “  THT  ’  y*  "  7^(1)  (7'?) 

The  associated  expansion  is  of  the  form 

♦(*>?;*)  -  X  +  vd5F  +  u5/3  <P(x*,y*)  +  •••  (7.8) 
where 

vd  «  vertical  flow  at  (x«xd,y»0)  - 

“  *0_(V°+)  • 

y 

The  equation  and  surface  boundary  condition  of 
tangent  flow  for  the  local  region  is 

(y+i)  VW  -  Vy*  ‘  °  (7  *9) 

V(x*’0+)  “  **  Fu(xd)  (7*10) 

Only  the  local  curvature  is  important  here.  A 
solution  must  be  found  with  a  shock  wave  and  this 
must  match  asymptotically  as  (x*-«,y*-») 

(x-x  J  -  0+)  to  the  outer  expansion.  ThiB 
matching  determines  u(e)  as  well  as  providing 
the  boundary  condition  for  <p. 
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